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ABSTRACT 

Recent  military  conflicts  in  Iraq  and  Afghanistan  have 
highlighted  the  wartime  effect  of  traumatic  brain  in¬ 
jury  (TBI).  While  it  is  not  clear  why  TBI  has  been  so 
prominent  in  these  particular  conflicts,  one  reason  may 
be  that  improvements  in  body  armor  have  led  to  in¬ 
creased  survivability  of  blasts.  Closed  traumatic  brain 
injury  covers  a  spectrum  of  central  nervous  system 
(CNS)  injuries  and  mechanisms,  but  it  is  broadly  char¬ 
acterized  as  mild  (mTBI),  moderate,  or  severe  TBI, 
with  mTBI  occurring  most  frequently.  Blunt,  ballis¬ 
tic,  and  blast  effects  may  all  contribute  to  CNS  injury, 
but  blast  in  particular  has  been  suggested  as  a  primary 
cause  of  military  TBI.  Little  is  currently  known  about 
the  effects  of  blasts  on  the  CNS;  injury  thresholds  have 
not  been  established,  and  even  direct  transmission  of 
the  non-linear  shock  wave  into  the  intra-cranial  cav¬ 
ity  and  brain  is  disputed.  In  this  study,  we  demon¬ 
strate  the  potential  for  a  blast  shockwave  to  directly 
affect  the  CNS  using  coupled  computational  fluid-solid 
dynamics  simulation.  The  model  includes  a  complex 
finite  element  model  of  the  head  and  intracranial  con¬ 
tents.  The  effects  of  threshold  and  lethal  blast  lung 
injury  were  compared  with  concussive  impact  injury 
using  the  full  head  model. 

1.  Introduction 

Current  military  operations  in  Iraq  (Operation  Iraqi 
Freedom)  and  Afghanistan  (Operation  Enduring  Free¬ 
dom)  have  brought  into  sharp  focus  military-related 
traumatic  brain  injury  (TBI).  TBI  is  a  significant 
cause  of  death  and  morbidity  in  the  0  to  40-year-old 
range  with  a  huge  direct  and  indirect  economic  im¬ 
pact  through  the  burden  of  care  imposed  on  family 
members  and  society  at  large  (Brazarian  et  ah,  2005; 
Bruns  and  Hauser,  2003).  Although  blasts  have  been 
suggested  as  a  primary  cause  of  military  TBI,  little 
is  currently  known  about  the  effects  of  blasts  on  the 
CNS.  Blast  effects  on  some  biological  tissues,  such  as 
the  lung,  are  documented  in  terms  of  injury  thresh¬ 


olds,  but  this  is  not  the  case  for  the  CNS.  Even  di¬ 
rect  transmission  of  the  non-linear  shock  wave  into 
the  intra-cranial  cavity  and  brain  is  disputed,  with  one 
suggested  mechanism  of  potential  CNS  injury  being  in¬ 
direct  transmission  of  the  blast  compression  wave  via 
the  cerebral  vasculature.  There  is  no  current  experi¬ 
mental  evidence  that  blast  waves  per  se  alter  the  CNS 
injury  cascade  of  TBI  compared  to  blunt  injury. 

A  blast  wave  is  a  pressure  wave  of  finite  amplitude 
resulting  from  an  atmospheric  explosion  that  releases 
a  large  amount  of  energy  in  a  short  period  of  time. 
In  addition  to  the  blast  wave,  the  chemical  energy  of 
the  explosion  is  released  into  thermal,  electromagnetic, 
and  possibly  into  kinetic  energy  of  solid  fragments  and 
surrounding  material  (Strehlow  and  Baker,  1976).  The 
non-linear  shock  wave  represents  a  large  discontinuous 
increase  in  pressure  and  density  in  the  air  which  prop¬ 
agates  radially  from  the  source  of  the  explosion  with 
the  potential  to  cause  large  loading  transients.  Impor¬ 
tantly,  an  explosive  event  is  stochastic  with  the  charac¬ 
teristics  of  a  three-dimensional  complex  flow  field  that 
may  be  altered  by  ambient  conditions  and  environ¬ 
mental  boundaries  that  can  result  in  wave  reflections 
and  up  to  eight-fold  field  intensification  (Cullis,  2001; 
Kambouchev  et  ah,  2007).  A  free  field  blast  wave  can 
be  most  simply  described  by  the  idealized  Friedlander 
waveform  with  a  rapid  rise  to  the  peak  pressure,  ex¬ 
ponential  fall-off  of  the  overpressure,  together  with  a 
relatively  prolonged  underpressure,  resulting  in  a  com¬ 
bination  of  compressive  and  tensile  components  as  it 
propagates  through  a  material  (Taylor,  1950;  von  Neu¬ 
mann,  1943).  The  interaction  of  a  blast  wave  with  a 
material  either  biological  or  non-biological  results  in 
the  development  of  stress  waves  which  have  longitu¬ 
dinal,  shear  (transverse),  and  shear  (Rayleigh)  wave 
components  (Kolsky,  1963).  Development  of  a  shear 
wave  where  the  direction  of  particle  motion  is  orthog¬ 
onal  to  the  direction  of  wave  propagation  results  in  ma¬ 
terial  shear  stresses  that  may  have  a  different  patho¬ 
logical  effect  on  the  anisotropic  structures  of  the  brain, 
especially  white  matter  tracts. 

Coupling  of  the  nonlinear  blast  shockwave  into 
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biological  tissue  results  in  significant  energy  depo¬ 
sition  at  high  strain  rates  in  fractions  of  millisec¬ 
onds.  The  relative  amount  of  energy  deposition  will  be 
most  simply  dependent  on  the  distance  from  the  blast 
source  (inverse  cube  law)  together  with  the  high  strain 
rate  tissue  material  properties.  Recent  data  suggests 
that  tissue  shear  response  may  exhibit  differing  states 
across  the  range  of  high  strain  rates  as  determined 
by  Kolsky  bar  experiments.  Similar  experiments  in¬ 
dicated  that  the  bulk  modulus,  while  non-linearly  re¬ 
lated  to  strain  rate,  did  not  exhibit  more  than  one 
state  (Saraf  et  ah,  2007). 

The  propagation  of  stress  waves  in  the  brain  and 
their  relation  to  TBI  have  been  studied  via  simulation 
under  impact  loading  (e.g.,  Zhang  et  ah,  2004;  Will- 
inger  and  Baumgartner,  2003).  However,  due  to  the 
characteristics  of  blast  waves,  it  is  expected  that  the 
loading  transients  under  blast  conditions  may  involve 
much  shorter  chact eristic  time  scales,  leading  to  higher 
strain  rates  than  seen  under  impact  conditions.  This 
may  set  in  motion  deformation  mechanisms  of  tissue 
response  that  are  not  operative  under  impact  condi¬ 
tions,  which  may,  in  turn,  lead  to  a  differing  spectrum 
of  blast-related  brain  injury. 

The  purpose  of  this  paper  is  to  develop  a  full  head 
model  and  simulation  of  the  fluid-solid  interaction  with 
a  blast  shockwave  and  to  interpret  these  events  with 
reference  to  the  Bowen  curves,  which  give  the  esti¬ 
mated  tolerance  to  a  single  blast  at  sea  level  for  a  70- 
kg  human  oriented  perpendicular  to  the  blast  (Bowen 
et  ah,  1968).  Estimates  from  the  Bowen  curves  suggest 
that  the  threshold  for  unarmored  parenchymal  lung  in¬ 
jury  at  standoff  of  0.6  m  is  80  psi  or  5.45  atm,  while 
400  psi  or  27.21  atm  results  in  1%  survival  from  lung 
injury.  250  psi  or  17.01  atm  is  the  LD$ 0  (lethal  dose, 
50%)  with  approximately  50%  survival  from  lung  in¬ 
jury  (Bowen  et  ah,  1968). 

In  order  to  examine  and  compare  the  effects  of 
blast  shock  waves  on  the  human  head,  simulations 
were  run  with  a  full  head  mesh  in  three  different  con¬ 
texts:  (1)  a  blast  with  overpressure  of  5.2  atm  or 
threshold  lung  injury,  equivalent  to  a  free  air  explo¬ 
sion  of  0.0648  kg  TNT  at  a  0.6  m  standoff  distance; (2) 
a  blast  with  overpressure  of  18.6  atm  or  near  the  LD$o 
of  lung  injury  survival,  equivalent  to  a  free  air  explo¬ 
sion  of  0.324  kg  TNT  at  a  0.6  m  standoff  distance;  and 
(3)  an  impact  between  a  head  traveling  at  5  m/s  and 
a  stationary,  immovable  boundary.  The  overpressures 
for  the  two  blast  simulations  were  selected  based  on 
the  Bowen  curves(Bowen  et  ah,  1968).  In  the  impact 
simulation,  it  was  determined  that  the  impact  velocity 
should  be  5  m/s  in  order  to  result  in  concussive  injury, 
based  on  comparable  impact  studies  in  the  literature 
(Casson  et  ah,  2008;  Zhang  et  ah,  2004). 


2.  Methods 

High  resolution  T1  MR  images  were  downloaded  from 
the  Montreal  Neurological  Institute  at  an  isotropic 
voxel  dimension  of  1  x  1  x  1  mm  (Collins  et  ah,  1998). 
These  images  were  merged  with  a  bone  windowed  CT 
of  the  head  allowing  skull  reconstruction  using  a  mu¬ 
tual  information  algorithm.  The  resulting  volume  set 
of  images  was  then  segmented  by  hand  into  topologi¬ 
cal  closed  regions  of  interest  using  Amira  and  then  ex¬ 
ported  as  VRML  files.  These  files  were  imported  into 
ICEMCFD,  and  unstructured  finite  element  meshes 
were  formed  using  Octree  and  Delaunay  tetrahedral 
mesh  generation  algorithms.  The  meshes  were  further 
refined  by  isolation  of  poorly  meshed  areas  and  poorly 
shaped  tetrahedra  prior  to  running  the  computational 
fluid-solid  dynamics  code  (Deiterding  et  ah,  2006).  It 
was  found  that  meshes  with  fewer  than  700,000  ele¬ 
ments  were  too  coarse  to  describe  the  intricate  geome¬ 
try  and  topology  of  some  anatomical  structures  of  the 
human  head  deemed  relevant  for  blast  injury  analysis. 
Previous  computational  models  have  been  reported  in 
the  literature  having  a  maximum  of  314,500  elements 
(Zhang  et  ah,  2001). 

In  order  to  balance  the  mesh  resolution  and  com¬ 
putational  requirements,  a  mesh  with  808,766  elements 
was  used  in  the  simulations.  Meshes  of  higher  defini¬ 
tion  will  be  evaluated  in  future  work  in  combination 
with  higher-fidelity  constitutive  models  and  proper¬ 
ties  of  tissue  reponse.  The  computational  model  dif¬ 
ferentiates  11  distinct  structures  characterized  by  me¬ 
chanical  function  rather  than  physiology:  ventricular 
cerebrospinal  fluid,  peri- ventricular  glia,  white  matter, 
gray  matter,  eyes,  venous  sinuses,  subarachnoid  cere¬ 
brospinal  fluid,  air  sinuses,  muscle,  skin  and  fat,  and 
diploic  skull  bone. 

3.  Material  Models 

The  constitutive  response  of  brain  tissue  encompasses 
a  variety  of  complex  mechanisms  including  nonlinear 
viscoelasticity,  anisotropy,  and  a  strong  rate  depen¬ 
dence  (Velardi  et  ah,  2006;  Shen  et  ah,  2006).  Sev¬ 
eral  investigations  have  focused  on  characterizing  this 
response  experimentally  and  on  developing  a  variety 
of  constitutive  models  capturing  this  behavior  (e.g., 
Miller  and  Chinzei,  2002;  Prange  and  Margulies,  2002). 
Owing  to  the  complexities  and  inherent  variability  as¬ 
sociated  with  biological  tissue,  there  is  usually  signif¬ 
icant  uncertainty  in  quantifying  tissue  response.  In 
consideration  of  these  limitations,  computational  mod¬ 
els  have  usually  favored  simpler  (i.e.  elastic)  mod¬ 
els  with  few  parameters  that  can  be  quantified  with 
less  uncertainty,  instead  of  more  sophisticated  models 
with  many  parameters  that  are  harder  to  estimate.  In 
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the  impact  TBI  modeling  work,  isotropic  elastic  mod¬ 
els  have  been  used  for  the  volumetric  response  while 
linear  viscoelastic  effects  have  been  considered  in  the 
shear  response  via  a  time-dependent  shear  modulus 
evolving  from  the  instantaneous  to  the  long-term  value 
(Belingardi  et  ah,  2005;  Willinger  and  Baumgartner, 
2003;  Zhang  et  ah,  2001).  The  relaxation  times  in¬ 
volved  are  on  the  order  of  tens  to  hundreds  of  millisec¬ 
onds  or  higher.  It  is  therefore  reasonable  to  expect 
that  deferred  deformation  or  stress  relaxation  due  to 
viscoelastic  effects  play  a  secondary  role  under  blast 
loading,  where  the  characteristic  times  seldom,  if  ever, 
exceed  a  few  milliseconds. 

On  this  account,  and  as  a  first  approximation,  we 
have  adopted  a  simplified  constitutive  modeling  strat¬ 
egy  emphasizing  the  effects  pertinent  to  blast  condi¬ 
tions,  specifically  the  description  of  the  pressure  wave 
propagating  through  the  brain  via  a  suitable  equation 
of  state.  To  this  end,  the  volumetric  response  of  brain 
tissue  has  been  described  by  the  Tait  equation  of  state 
with  parameters  adjusted  to  fit  the  bulk  modulus  of 
the  various  tissue  types.  The  deviatoric  response  has 
been  described  via  a  large-deformation,  neo-Hookean 
elastic  model  with  properties  adjusted  to  fit  reported 
values  of  the  instantaneous  shear  modulus.  The  Mie- 
Gruneisen/Hugoniot  equation  of  state  was  used  to  de¬ 
scribe  the  volumetric  response  of  the  skull.  The  consti¬ 
tutive  properties  of  the  tissues  were  determined  from 
a  literature  review. 


3.1  Mie-Gruneisen/Hugoniot  Equa¬ 
tion  of  State 


The  shock  response  of  many  solid  materials  is  well  de¬ 
scribed  by  the  Hugoniot  relation  between  the  shock 
velocity  Us  and  the  material  velocity  Up  of  the  sim¬ 
ple  form  (D.S.  Drumheller,  1998;  M.A.  Meyers,  1994; 
Zel’dovich  and  Raizer,  1967): 

Ua  =  Co  +  sUp  (1) 


In  this  expression,  Co  and  s  are  material  parameters 
that  can  be  obtained  from  experiments.  Here  the  val¬ 
ues  for  Co  and  s  are  the  same  as  those  used  for  the 
skull  in  Taylor  and  Ford  (2006).  By  considering  Equa¬ 
tion  (1)  and  conservation  of  mass  and  momentum  in 
a  control  volume  at  the  shock  front,  the  final  pressure 
can  be  calculated  explicitly  as  a  function  of  the  Ja¬ 
cobian  behind  the  shock  front  Jh  and  the  reference 
density  ahead  of  the  shock  p0  (D.S.  Drumheller,  1998; 
M.A.  Meyers,  1994;  Zel’dovich  and  Raizer,  1967): 


PoClO-  JH) 


(2) 


where  Jh  is  related  to  the  density  p#,  the  specific 
volume  Vh,  or  the  deformation  gradient  tensor  F#, 


defined  behind  the  shock  front,  by: 

Jh  =  —  =  rry-  =  det(F H)  (3) 

Ph  Vo 

The  relation  Equation  (2),  called  the  shock  Hugoniot, 
relates  any  final  state  of  density  to  its  corresponding 
pressure.  The  deformation  path  taken  by  the  mate¬ 
rial  between  the  initial  state  (Pq,  Vo)  and  the  final 
state  (Ph,  Vh)  is  then  defined  by  the  Rayleigh  line:  a 
straight  line  in  the  (<ti,  V)  plot,  where  G\  is  the  axial 
stress  in  the  shock  direction  (D.S.  Drumheller,  1998; 
M.A.  Meyers,  1994).  The  parameters  used  in  the  sim¬ 
ulations  are  given  in  Table  1. 

Structure  p(kg/m 3)  E(Pa)  v  Co  S 

Skull  1412  6.5e9  022  1.85e3  094 

Table  1:  Material  Parameters  for  the  Skull 

3.2  Tait  Equation  of  State 

The  Tait  equation  of  state,  which  is  commonly  used  to 
model  fluids,  is  given  by  (Thompson,  1972): 

P  =  B[(p/p0)r°+1  ~  1]  (4) 

where  B  and  To  are  constants.  The  Tait  equation  of 
state  was  selected  to  describe  the  volumetric  response 
of  all  head  structures  except  the  skull  because  the  re¬ 
maining  structures  largely  consist  of  water.  To  ob¬ 
tain  the  necessary  parameters,  T0  was  taken  to  be  the 
value  for  water,  6.15,  appropriate  bulk  modulus  val¬ 
ues  K  were  selected  from  the  literature,  and  B  was 
computed  for  each  structure  using  the  relation: 

K0  =  B  x  (T0  +  1)  (5) 


3.3  Deviatoric  Elasticity 

In  addition,  all  of  the  structures  were  assumed  to  obey 
an  elastic  deviatoric  behavior.  To  this  end,  the  equi¬ 
librium  stress-strain  relation  was  extended  to  include 
the  elasticity  of  the  material  in  shear,  with  the  result 
(Cuitino  and  Ortiz,  1992;  Ortiz  and  Stainier,  1999): 

crB  =  -PI  +  J-1Fe[p(logVCP)dev\FeT  (6) 

where  Ce=FeT Fe  is  the  elastic  Cauchy-Green  defor¬ 
mation  tensor,  logVCP  is  the  logarithmic  elastic  strain, 
p  is  the  shear  modulus,  and  the  pressure  P  follows 
from  the  equations  of  state  defined  in  Equations  (2) 
and  (4).  Explicit  formulae  for  the  calculation  of  the 
exponential  and  logarithmic  mappings,  and  the  calcu¬ 
lation  of  their  first  and  second  linearizations,  has  been 
given  by  Ortiz  et  al.  (2001).  For  the  simulations,  val¬ 
ues  of  p  were  selected  from  the  literature.  These  values 
are  listed  in  Tables  1  and  2. 
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Structure 

p(kg/mn3) 

K  (Pa) 

B 

r0 

g(Pa) 

Ventricle 

1040 

2.19e9 

3.063e8 

6.15 

2.25e4 

Glia 

1040 

2.19e9 

3.063e8 

6.15 

2.253e4 

White  Matter 

1040 

2.19e9 

3.063e8 

6.15 

2.253e4 

Gray  Matter 

1040 

2.19e9 

3.063e8 

6.15 

2.253e4 

Eyes 

1040 

2.19e9 

3.063e8 

6.15 

2.253e4 

Venous  Sinus 

1040 

2.19e9 

3.063e8 

6.15 

2.013e4 

CSF 

1040 

2.19e9 

3.063e8 

6.15 

2.253e4 

Air  Sinus 

1040 

2.19e9 

3.063e8 

6.15 

2.253e4 

Muscle 

1100 

3.33e6 

4.662e5 

6.15 

3.793e4 

Skin/Fat 

1040 

3.47e7 

4.866e6 

6.15 

5.880e6 

Table  2:  Material  Parameters  for  Head  Structures 


4.  Simulation  Details 

The  two  blast  fluid-solid  interaction  simulations  were 
run  on  14  processors  for  the  solid  and  6  processors  for 
the  fluid,  with  two  levels  of  subdivision  allowed  for  the 
fluid.  The  lower  region  of  the  head  was  fixed  in  order 
to  avoid  the  blast  engulfing  the  bottom  of  the  head, 
where  the  neck  would  ordinarily  be  located.  The  solid- 
only  impact  simulation  was  run  on  20  processors. 

5.  Results 

5.1  5.2  atm  Simulation 

Figures  l(a)-(f)  illustrate  the  propagation  of  the  com¬ 
pressive  blast  wave  through  the  coronal  sections  of  the 
head  in  the  5.2  atm  simulation.  The  blast  wave  is  in¬ 
cident  on  the  right  temporal  region.  The  compressive 
wave  is  seen  propagating  through  the  cranial  cavity 
from  the  right  to  the  left  with  some  minor  reflection 
from  the  left  side  of  the  cavity,  leading  to  a  pocket  of 
concentrated  pressure  in  the  skull  on  the  right  side  of 
the  head. 

The  maximum  tensile  and  compressive  pressures 
were  extracted  and  plotted  for  each  time  step  for  each 
of  the  11  distinct  structures.  These  curves  are  given  in 
Figures  2 (a) -(f)  for  the  blast  and  impact  simulations. 
In  the  5.2  atm  simulation,  the  maximum  compressive 
pressure  reached  was  6.5  MPa  at  0.00045  s,  and  the 
maximum  tensile  pressure  reached  was  0.89  MPa  at 
0.00048  s.  It  is  evident  from  the  curves  that  the  var¬ 
ious  structures  experienced  different  time  courses  and 
different  pressure  intensities.  The  highest  compressive 
pressures  were  experienced  by  the  skull,  muscle,  and 
subarachnoid  CSF,  while  the  highest  tensile  pressures 
were  experienced  by  the  air  sinuses,  gray  matter,  sub¬ 
arachnoid  CSF,  skull,  and  white  matter.  In  this  sim¬ 
ulation,  the  nodes  that  experienced  the  highest  pres¬ 
sures  were  all  located  on  the  right  side  of  the  head, 
in  the  concentrated  pocket  of  stress  created  from  the 
reflection  of  the  blast  wave  from  the  left  side  of  the 


head. 

5.2  18.6  atm  Simulation 

At  the  LD$o  blast  wave  overpressure  of  18.6  atm,  the 
compressive  propagation  is  seen  in  Figure  3  with  gen¬ 
eration  of  differential  stresses  as  the  wave  propagates. 
Wave  reflection  is  more  apparent  at  this  overpressure 
as  seen  in  the  multiphasic  curves  of  Figures  2(c)  and 
(d).  The  simulation  reached  a  maximum  compressive 
pressure  of  39  MPa  at  0.00042  s  in  the  skull  and  mus¬ 
cle  and  a  maximum  tensile  pressure  of  4.5  MPa  at 
0.00041  s  in  the  skull  and  gray  matter.  High  compres¬ 
sive  pressures  were  experienced  by  the  skull,  muscle, 
subarachnoid  CSF,  gray  matter,  and  skin/fat  mate¬ 
rial  elements.  The  skull  experienced  peak  compres¬ 
sive  pressure  first,  followed  by  the  subarachnoid  CSF 
and  muscle,  and  then  gray  matter.  High  tensile  pres¬ 
sures  were  experienced  by  the  skull,  subarachnoid  csf, 
gray  matter,  and  white  matter.  Compared  to  the  5.2 
atm  simulation,  the  18.6  atm  simulation  experienced 
pressures  that  were  4-6  times  higher  than  those  expe¬ 
rienced  in  the  5.2  atm  simulation,  and  the  pressures 
and  stresses  peaked  earlier  in  the  18.6  atm  simulation. 
The  locations  in  the  head  that  experienced  the  highest 
stresses  were  very  similar  in  both  the  5.2  atm  and  18.6 
atm  simulations. 

5.3  Impact  Simulation 

The  5  m/s  impact  simulation  ran  to  .000634  s,  reach¬ 
ing  a  maximum  compressive  pressure  of  27.2  MPa  and 
a  maximum  tensile  pressure  of  7.1  MPa.  The  impact 
was  delivered  in  the  mid-coronal  plane  in  a  lateral  di¬ 
rection,  from  left  to  right.  Significant  differences  are 
seen  between  blast  and  impact  injury,  specifically  in 
the  monotonic  form  of  the  impact  compression  and 
tension  curves. 

High  compressive  pressures  were  experienced  by 
the  skull,  subarachnoid  CSF,  muscle,  skin/fat,  and 
gray  matter,  with  the  highest  compressive  pressure  be¬ 
ing  experienced  by  the  skull.  High  tensile  pressures 
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Figure  1:  Propagation  of  blast  compressive  wave  through  cranial  cavity  in  5.2  atm  simulation,  viewed  through  mid¬ 
coronal  plane  at  0.132-0.685  ms.  The  propagation  of  the  blast  wave  is  from  left  to  right  in  the  external  observer 
frame  of  reference.  Scale  is  from  0  to  500  kPa. 


were  experienced  by  the  skull,  muscle,  skin/fat,  and 
gray  matter,  again  with  the  highest  tensile  pressure 
being  experienced  by  the  skull.  The  magnitudes  of 
the  maximum  pressures  reached  in  the  impact  sim¬ 
ulation  are  comparable  to  the  maximum  pressures 
reached  in  the  18.6  atm  simulation,  suggesting  that 
blasts  that  would  result  in  50%  lethality  due  to  unar¬ 
mored  parenchymal  lung  injury  would  lead  to  concus- 
sive  injury  in  the  brain. 

CONCLUSIONS 

The  use  of  personal  protective  equipment  (PPE),  par¬ 
ticularly  body  armor,  is  considered  to  be  protective 
against  lethal  blast  and  fragmentation  lung  injury. 
The  goal  of  this  work  was  to  use  the  standardized 
Bowen  curves  for  threshold  and  50%  lethal  lung  in¬ 
jury  to  examine  the  effects  of  similar  exposures  on 
the  head  and  to  compare  them  with  impact  decelera¬ 
tions  such  as  those  often  seen  within  sports  concussions 
(Zhang  et  ah,  2004).  The  analysis  indicates  that  injury 
from  a  blast  resulting  in  50%  lethal  unarmored  lung 
injury  is  associated  with  impact  injury  that  may  result 
in  concussion.  The  American  Academy  of  Neurology 
grades  concussions  from  1-3,  with  grade  3  implying 
a  loss  of  consciousness  and  grades  1  and  2  implying 
transient  confusion  (American  Academy  of  Neurology, 
1997).  In  the  military  context  such  transient  confusion 
may  lead  to  a  loss  of  situational  awareness  within  the 
combat  context,  with  potentially  devastating  results. 
These  results  need  to  be  experimentally  validated  and 
are  best  considered  by  regarding  TBI  to  occur  across 
a  strain-rate  continuum  with  impact  injury  occurring 
at  lower  strain  rates,  ballistic  injury  at  intermediate 
to  higher  strain  rates,  and  blast  injury  at  the  highest 
strain  rates  (Moore  et  ah,  2008).  Such  validation  re¬ 
quires  both  experimental  measurement  but  also  appro¬ 
priate  models  of  tissue  constitutive  properties  across 
the  strain-rate  continuum.  The  constitutive  models 
used  in  the  current  simulations  -  the  Tait  and  Mie- 
Grueneisen/Hugoniot  equations  of  state  -  while  likely 
to  be  refined  by  further  experimental  work,  represent 


good  first-order  approximations  that  are  likely  to  hold 
at  impact  and  blast  strain  rates. 

The  three  contexts  compared  in  this  paper  were: 
(a)  a  peak  blast  overpressure  of  5.2  atm  corresponding 
to  the  Bowen  threshold  for  unarmored  lung  parenchy¬ 
mal  injury;  (b)  a  peak  blast  overpressure  of  18.6  atm 
corresponding  to  the  Bowen  LD^q  for  unarmored  lung 
parenchymal  injury;  and  (c)  impact  deceleration  of  the 
full  head  model  from  a  velocity  of  5  m/s  to  0  m/s  fol¬ 
lowing  impact  with  a  stationary  immovable  boundary. 
These  conditions  were  chosen  specifically  to  consider 
whether  there  may  be  calculable  evidence  for  concus- 
sive  brain  injury  following  primary  blast  exposure  at 
a  specific  stand-off.  This  would  appear  to  be  the  case 
under  the  strict  model  conditions  applied  and  suggest 
at  a  minimum  that  such  a  potential  threat  deserves 
considered  further  investigation. 

Currently,  head  PPE  tends  to  be  optimized  for 
impact  or  ballistic  protection,  with  little  or  no  con¬ 
sideration  of  blast  mitigation  or  protection.  The  engi¬ 
neering  optimization  issues  required  to  maximum  head 
PPE  should  not  be  underestimated,  but  consideration 
of  optimizing  two  domains  such  as  impact  and  blast 
injury  may  be  obtainable. 
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